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We present a new method based on the scattering technique to investigate fluctuation-induced 
forces at a fluid interface. The scattering approach, well suited to the study of many body systems 
of arbitrary geometries, is augmented to include boundary fluctuations. Using this method, we 
study the deviation of the total fluctuation-induced interaction from the sum of pairwise energies 
for three colloidal particles. We consider both frozen and fluctuating colloids and obtain a very 
good agreement between analytical and numerical results. We find a marked difference in the three 
body fluctuation-induced free energy between the frozen and fluctuating colloids, both in sign and 
relative size. 



With the construction of micro- and nano-eletro- 
mechanical devices (MEMS and NEMS) [TH3] and the 
desire to understand biological processes in greater de- 
tails, it has become necessary to have a complete under- 
standing of fluctuation- induced forces. These forces were 
first predicted by Casimir in 1948 who showed that con- 
finement of zero point energy fluctuations of the electro- 
magnetic field between two conducting plates in vacuum 
gives rise to long ranged attractive interactions [H [5]- 
Thirty years later, Fisher and de Gennes noted that im- 
posing constraints on the thermal, rather than quantum, 
fluctuations in fluids in the vicinity of their critical points 
also results in similar long-ranged forces [B]. Since then, 
the thermal versions of Casimir forces have been inves- 
tigated in various theoretical and experimental studies 

[THIS]- 

In particular, due to the industrial applications and im- 
portance in biological systems, the thermal fluctuation- 
induced forces between colloidal particles have been stud- 
ied within the last few years. For example, the experi- 
ments of Hertlein et. al. reveal the presence of repulsive 
and attractive fluctuation-induced forces on colloidal par- 
ticles immersed in a near-critical binary mixture in the 
vicinity of a wall [11] . The thermal fluctuation- induced 
force also plays an important role in the effective inter- 
action between colloidal particles trapped between two 
fluid phases [2pJ[2T]- 

There have been previous works using various tech- 
niques to study fluctuation-induced forces at interfaces 
with interesting results [T21 - H9] . There are, however, a 
number of complications associated with each of these 
methods. The most common is the difficulty of calculat- 
ing the interaction between colloids at very small separa- 
tions and matching it with the Derjaguin approximation 
[22] . Moreover, extending the previous methods to vari- 
ous geometries and many particle systems can be tedious. 

In this work we study the interaction between col- 
loids trapped at a fluid interface using the scattering 
formalism 23-25], which is widely used to obtain quan- 
tum fluctuation-induced forces. The scattering formalism 
presents a straightforward approach for treating various 



geometries and many body systems and is numerically 
faster compared to previous methods. 

Furthermore, the scattering method simplifies the 
fluctuation-induced calculations by writing the energy 
into the translation matrices (U) and scattering matri- 
ces or T-matrices (T). The translation and T-matrices 
describe how a fluctuation propagates through the field 
between the objects and how fluctuations interacts with 
the objects, respectively. 

The important difference between the QED Casimir ef- 
fect and the fluctuation-induced forces between colloids 
at an interface is that the colloids and their boundaries 
can fluctuate. This work extends the scattering method 
to include these effects, by separating out the interface 
and colloid Hamiltonians. This new extension can sim- 
plify the problem by separating the properties of the col- 
loid, as well as the type of fluctuations, from the calcula- 
tion of the T-matrix. Note, in contrast to previous works 
[T2TH61 [T9] , the calculation of the T- and the translation 
matrices are completely decoupled. 

In order to demonstrate the effectiveness of this new 
technique, we calculate the fluctuation-induced forces be- 
tween three spherical Janus colloids at an interface. We 
study three different fluctuation scenarios: (i) colloids 
frozen at the interface; (ii) bobbing colloids that fluctuate 
only vertically; and (iii) bobbing and tilting colloids that 
both fluctuate vertically and tilt side to side. Surpris- 
ingly we find different three body behaviors for fixed and 
fluctuating colloids. The three body effect for the fixed 
colloids is repulsive and comparable to the two body in- 
teraction, while for fluctuating colloids is attractive and 
comparable only at short separations. 

In the following, we briefly explain the derivation of 
the scattering formalism for colloids at fluid interfaces. 
Moreover, we show some numerical and analytical results 
obtained with this technique. Note that due to the rapid 
communications format, we only outline the important 
steps. The details of the derivations and calculations are 
left for a more complete exposition [26] . 

We consider an infinite interface between two fluid 
phases characterized by the surface tension a. At equi- 
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librium, the interface is flat and placed at z = 0. We sup- 
pose that the thermal height fluctuations of the interface 
are small without overhanges and bubbles. Therefore, we 
express the interface height profile in the Monge repre- 
sentation, i.e. z — w(x). The capillary wave Hamiltonian 
[T51 [TB] then becomes 
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where A c = \fa~f{Apg) is the capillary length with Ap = 
pi — p2, the difference in the mass density of the two fluid 
phases. 

We now introduce N colloidal particles at the inter- 
face between two fluid phases. The colloid Hamiltonian 
will be defined as the energy costs associated with the 
insertion of a colloid, 
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where hi is the height of the center of mass of the colloid 
i with respect to the equilibrium interface, and cr^ini) is 
the surface tension between the colloid i and the fluid 
phase I (II). The contact line field /i(x) is the height at 
which the fluid interface intersects the colloid, and is ex- 
tended into the interior of the colloid. The cross sections 
of the colloid i with the equilibrium and fluctuating in- 
terfaces are denoted by f2J cf and respectively, and 
A£l = £! rc f — f2 is the change in the projected area. The 
surface area of the colloid in contact with the fluid phase 
I(II) is represented by ^4i(n) where AAij is the change 
in the area of the colloid i in the fluid phase I. Note that 
for Janus particles the total area of the colloid does not 
change AAi — — A An. In addition, the gravitational en- 
ergy due to the colloid height hi is negligible and can be 
ignored in Eq. ^ because hi/\ c <gc 1. 

The first term in Eq. ^ is associated with excluding 
the projected areas fl; from the interface since the in- 
tegration domain in Eq. ([I]) covers the whole R 2 space. 
The second term represents the energy costs related to 
the change of the projected area for a colloid in the fluc- 
tuating vs reference interface, and the third term is the 
energy costs due to the change in the area of the colloids 
in the phase I and II. 

The total Hamiltonian of the system is the sum of the 
interface and N colloid Hamiltonians, Utot = "HintM + 
/ H l co i[fi, hi] , and the partition function Z is the sum 
over all interface and colloid configurations 
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with C representing the constraints imposed by the col- 
loidal particles on the interface height fluctuations. At 



the surface of the colloid i the interface height field u 
is constrained by the contact line field /,;. A Dirac 
delta functional is used to implement the constraints, 
S[u -/,•] = / n,=i Vipi exp[i J Qi d 2 x tp(u - fi)}, where we 
have introduced the auxiliary field ipi(x). Equation ^ 
can be separated as 



N 
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where Z- lnt and Z co \ correspond to the partition functions 
related to the interface and colloid fluctuations respec- 
tively. 

The integration over the scalar field u in the interface 
partition function can be performed and results in 
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where Co is a constant and 
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In Eq. ([6]), G(x, x') is the Green's function of the cap- 
illary wave equation for the free interface and can be 
expanded in terms of the solutions to the capillary wave 
equation, i.e. (—V 2 + A 2 )u(r) = 0. Using separation 
of variables in a coordinate system centered on a col- 
loid, the capillary wave equation has two solutions: the 
solution regular at origin, which corresponds to the inci- 
dent field in the scattering method 4> l ^ c (r); and a solution 
that dies off at infinity, which corresponds to the scat- 
tered field (r). The expansion of the Green's function 
G(x, x') and the auxiliary field ipi in terms of these solu- 
tions are G(x,x') = £ a c a <^ nc (x)<f Q ct (x') and ^(x) = 
J2 a <3/ a ^Q C (x) , respectively, with c Q the expansion coef- 
ficient and VPa the multipole coefficient. Plugging these 
equations into Eq. we find Gy = J2 a p c a~® ia^a^ sP 
for i ^ j and Gu = - J2 a fi ■ cp^ia^apl^Vip for i = j. 
The translation matrix U tJ describes the couplings be- 
tween the partial waves on distinct colloids. The scat- 
tering amplitude matrix TP relates the incident and scat- 
tered fields, and depends only on the shape of the colloid. 
Inserting the expressions for Gij in Eq. Q, we find 
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The matrix (C) in Eq. ^ is a diagonal matrix with the 
elements that are the coefficient of the Green's function 
expansion, and \&i is a vector with elements that are the 
multipole moments associated with the particle i. 
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We evaluate the colloid partition function following the 
same steps as the ones we took to calculate the interface 
partition function. The colloid Hamiltonian in Eq. ([2| 
is expanded to be quadratic in the contact line field /j, 
and then the field is expanded in terms of the incident 
solutions, 

^(x) = 5> Q C c (x). (8) 

a 

After expansion the colloid partition function is 

N r N 

Z co , = / [] VPi exp - — Y, ?1 H col P, 

J i=l IB j=1 

-^(pt^ + p^t)^ (9) 

where the exact form of H col depends on the type of col- 
loid fluctuations and will be presented below. We now 
insert Eqs. ^ and Q into Eq. Q to obtain the total 
partition function Z. The resulting partition function is 
a Gaussian integral over the multipoles and P. The 
fluctuation-induced free energy at temperature T can be 
then obtained through E/ksT — — In (Z/Zoo), with Z^ 
corresponding to the partition function with all the col- 
loids placed at infinite distance from each other. The 
fluctuation-induced energy as a function of the T-matrix 
and the translation matrix U then reads as 



JL = Ilndet(l-TU), 
k B ± 2 
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where U, 



and T, 



modified T-matrix of the colloid i given by 



TSij with T the 



T i =T i -T i [CT i C- 1 + -Cff co ,] *T 
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Note while T depends only on the shape of the colloid, T 
also depends on H co i which in turn depends on the spe- 
cific fluctuations that the colloid undergoes. Equations 
(10) and (11) can be used to obtain the fluctuation- 
induced interaction between N colloidal particles at a 
fluid interface. 



We now employ Eq. ( 10 ) to calculate the fluctuation- 



induced free energy between three identical spherical col- 
loids with radius R trapped at an interface. For simplic- 
ity we assume that the colloids are Janus particles, i.e. 
the contact line is pinned to the colloid surface such that 
the area of the colloidal particles in two fluid phases does 
not change, i.e. AAij — 0. In addition, we set the equi- 
librium contact angle to tt/2. Using the solutions to the 
Hclmholtz equation in polar coordinates, the T-matrix 
for colloidal particles with circular cross section at the 
flat interface reads 



TV 
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K m (R/\ c ) 
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where I m and K m are the modified Bessel functions of 
the first kind and second kind, respectively. In order to 



obtain the fluctuation-induced energy given by Eq. (10), 
it is necessary to calculat e H ^i or more specifically the 

To find H ro i wc assume 
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modified T given in Eq. 
that the colloids are: (i) frozen, allowing no fluctuation; 
(ii) bobbing, allowing the colloids to fluctuate up and 
down; and (iii) bobbing and tilting, allowing the colloids 
to fluctuate both up and down as well as tilt side to side. 

For a frozen colloid at the interface, the colloidal par- 
ticles are exposed to an infinite potential barrier to fluc- 
tuate, therefore the elements of the matrix Tl co \^ tend 
to infinity and from Eq. (JTIJ T = T. For bobbing col- 
loidal particles, the projected area does not change and 
therefore Ail = in Eq. Substituting Eq. Q into 
Eq. (§, we can calculate the matrix H co i. In the limit 
R/X c -C 1, we find H co i ps —ikjR 2 /\ c for m — and 
H co i pa oo for \m\ > 0. From Eq. (11) wc immediately 
find that for m = 0, To ~ and for m 7^ 0, T m pa T m . 

For colloids that both bob and tilt, the projected 
area changes with fluctuations, i.e. Ail ^ 0. After 
calculating the change in the projected area for small 
tilt angles, using Eq. Q in the limit R/X c <C 1 we 
find H col — ■naR 2 j\ 2 cl 0, oo for m = 0, |m| = 1 and 
|m| > 1 respectively. Therefore, we have T ,T±i ss 0, 
for \m\ =0,1 and T TO pa T m otherwise. 

To calculate the fluctuation-induced energy given in 



Eq. ( 10 ), we also need to calculate the translation matrix 



in polar coordinate, which can be obtained by using the 
Graff's addition theorem for Bessel functions [57]. We 



find U 



y _ 



(-1)* 



t(m—m 



>(d/\ ) with 4>i 



the angle between the coordinate systems fixed to the 
colloids i and j. 

Employing the translation and T-matrices described 
above, we write 



i<j=i 
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where the first term is sum over the two-body interactions 
and the second term £"123 is the free energy due to the 
three body effect. For three fixed colloids at the interface, 
the asymptotic three body effect dij/Ri 3> 1 in the limit 
A c 3> Ri , dij reads 
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colloid i and dij the center-to-center distance between 



the colloids i and j. Note that the energy in Eq. (14) 
depends only on the distance between the colloids. 
For three bobbing colloids at the interface, the asymp- 
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FIG. 1. (Color online) Three body effect for colloids on a line 
with di2 = di3 = d.23/2 and on an equilateral triangle with 
di2 = di3 = d,23: The ratio of the three body energy £123 to 
the absolute value of the two body energy £23 vs. the surface- 
to-surface separation of the colloids 1 and 2, h/R. The dashed 
and dott ed lines represent the asymptotic results of Eqs. ( 15 I 
and ( 16 1 for the colloids in the linear and triangular configu- 



rations, respectively. Blue squares (linear configuration) and 
diamonds (triangular configuration) denote the numerical re- 
sults for the bobbing and tilting colloids. Red circles (linear 
configuration) and triangles (triangular configuration) show 
the numerical results for the bobbing only colloids. Inset: 
Three body effect for fixed colloids in the linear and trian- 
gular configurations. The dashed and dotted lines show the 



asymptotic result in Eq. (14). The black circles (linear config- 
uration) and red squares (triangular configuration) represent 
the numerical results for the fixed colloids. 
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Equation ( 15 1 is in agreement with the result in Refs. (TTJ 
ITS] and shows that the three body asymptotic energy is 
much smaller than the two body interaction d~ A . There- 
fore, at large separations, for bobbing colloids the sum of 
pairwise energies gives a very good approximation to the 
fluctuation-induced energy. 

For three bobbing and tilting colloids, the three body 
energy in the limit A c 3> Ri , dij reads 
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In this case the two body interaction scales with d~ 8 , and 
as such the pairwise energy summation again provides 
a good approximation for the total fluctuation-induced 
energy. 

Figure [T] depicts the ratio of the three body effects 
to the absolute value of the two body energies vs. the 



surface-to-surface separation, h/R — d^/R— 2, for two 
different configurations: colloids sitting on a line {dvz = 
d%3 = c?23/2) or at the vertices of an equilateral triangle 
(d±2 = d\3 = ^23)- As shown in the figure, there is 
a very good agreement between the asymptotic energies 



(dashed and dotted lines) given by Eqs. ( 14 1-( 16 ) and the 
numerical results (symbols) calculated by Eq. (10). 

As illustrated in the figure, the three body effect £123 is 
very small at large separations for bobbing only (red cir- 
cles and triangles) and bobbing and tilting colloids (blue 
squares and diamonds) in both configurations. In con- 
trast, the inset in Fig. [T] shows that the three body ef- 
fect energy £123 for fixed colloids is comparable to their 
two body interaction with the opposite sign, i.e., there 
is a strong deviation from the pairwise summation for 
the case of fixed colloids in both linear and triangular 
configurations, see the circles and squares in the inset. 
Quite interestingly, since the three body effect for fixed 
colloids is positive (repulsive) , the total energy is smaller 
than the sum of pairwise added interaction energies, i.e. 

£tot < £\2 + £13 + £23- 

In summary, we have presented an effective method 
to study the fluctuation-induced forces at fluid interfaces 
using a modified scattering method, in which the separa- 
tion of the interface and colloid free energies enables us to 
treat various interface and colloid fluctuations. The new 
approach has many advantages, such as the treatment of 
many colloid systems, quick numerical calculations, and 
the treatment of different colloid fluctuations. 

As an example, we studied the three body interaction 
between frozen and fluctuating spherical Janus colloids 
and found a very interesting difference between them. 
For both frozen and fluctuating colloids the two body 
interactions are attractive. However, for frozen colloids 
the three body interactions are repulsive weakening the 
interaction compared to the pairwise summation, while 
for fluctuating colloids the three body interactions are 
attractive strengthening the interaction. In addition, we 
found that the three body interaction for the frozen col- 
loids is non-negligible for the complete range of interac- 
tions, while the three body interactions for the fluctuat- 
ing colloids are negligible for larger separations. 

Due to the non-trivial behavior of these forces, a better 
knowledge of them will shed light on structure formation 
and crystallization phenomena at the fluid interfaces. 
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